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We show that excitons and free carriers in K and K’ valleys of transition metal dichalcogenide 
monolayers can be entangled with respect to their valley degree of freedom by absorbing linearly 
polarized single photons. This effect does not require any interaction between K and K’ excitons 
in contrast to conventional mechanisms of entanglement that are mediated by coupling between 
quantum systems (e.g. entanglement of photons in nonlinear optical interactions). The valley en¬ 
tanglement of excitons and free carriers can be verified by measuring the polarization of their 
photoluminescence or fluctuations of the photocurrent under an applied in-plane DC bias. 

PACS numbers: 


I. INTRODUCTION 

Transition metal dichalcogenides (TMDCs) of the composition MX 2 , where M = Mo or W and X = S, Se, or Te, 
are well known materials that have recently seen the resurgence of interest. This interest was largely driven by two 
discoveries. First, although bulk M 0 S 2 is an indirect-gap material, the monolayer M 0 S 2 turns out to have the direct 
band gap located in energy-degenerate K and K' points at the corners of the hexagonal Brillouin zone. This results 
in a dramatic increase of the band-edge photoluminescence (PL) yield by more that a factor of 10"^ [T]. Second, the 
combination of spatial inversion symmetry breaking with strong spin-orbit coupling leads to a valley-contrasting spin 
splitting of the valence-band edge, which gives rise to valley-dependent optical selection rules for interband transitions 
[2]. Namely, free carriers and excitons in the K and K' valleys are coupled to photons of the same energy but opposite 
helicity: left-hand circular (L) polarization for the K point and right-hand circular (R) polarization for the K' point. 
When excited by a circularly polarized light, M 0 S 2 monolayer shows PL with the same polarization as the excitation 
light, indicating that the valley polarization of excitons is preserved longer than the recombination time, at least at 
low temperatures [3H5]. When excited by a linearly polarized light, WSe 2 shows a high degree of linear polarization 
in the PL of neutral excitons, indicating that the inter-valley phase coherence survives in the processes of exciton 
formation and recombination [6]. Under an applied bias, valley-polarized TMDC monolayers exhibit valley and spin 
Hall effects [2] which enables optoelectronic devices based on the valley degree of freedom, e.g. valley Hall effect 
transistors [3- 

These recent results are very exciting as they suggest that the valley index degree of freedom in TMDC monolayers 
can serve as a robust information carrier. Excitons or free carriers in different valleys can be selectively manipulated 
by radiation in the convenient visible frequency range around 2 eV and by the in-plane DC electric field. A tanta¬ 
lizing question with implications for quantum information is whether and how one can achieve quantum mechanical 
entanglement of free carriers or excitons with respect to the valley index. 

The entanglement of two quantum systems or ensembles is usually generated as a result of coupling between them. 
This coupling can be mediated by classical electromagnetic fields; see e.g. [HHin]. However, a classical field cannot 
entangle non-interacting systems unless it directly couples them through a nonlinear optical process such as parametric 
frequency conversion or four-wave mixing. For example, the experiment in which employed a classical field cannot 
lead to entanglement of excitons (see Appendix A). 

At the same time, one can entangle non-interacting quantum systems by coupling them to a quantum field. Here we 
consider the optical excitation of electron-hole or neutral exciton states near the band gap of a MX 2 monolayer in two 
valleys K and K’ with opposite valley indices. We will consider entanglement of excitons for definiteness, although 
our analysis below works equally well for free electron-hole pairs. The binding energy of excitons in MX 2 monolayers 
is as high as several hundred meV, and the PL signal is dominated by excitons even for above band gap excitation. 
We will label two different valley states by an up or down pseudospin direction, assuming that (f) and (1) excitons 
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can be excited only by an R- or L-polarized field, respectively: 


Ei^,L oc e± 


x±iy 





FIG. 1: (Color online) A sketch of excitons in K and K’ valleys of MX 2 monolayers interacting with linearly polarized single 
photons. To ensure the interaction with one photon at a time, the power of incident radiation should not exceed hujAiu, where 
Auj is the radiation bandwidth. 


The main result of the paper is that absorption of linearly polarized single photons by a MX 2 monolayer in a cavity 
(Sec. II) or from a flux of incident photons (Sec. Ill) gives rise to an efficient entanglement of up and down excitons, 
i.e. of the valley degree of freedom; see Fig. 1. The intuitive physical explanation of this effect is that excitons in 
the K and K’ valleys of a MX 2 monolayer interact with a linearly polarized single photon as if it were an entangled 
R/L photon state. In fact, one can rigorously prove (see Appendix B) that a single photon state which is a factorized 
product state in a linearly polarized basis is equivalent to an entangled state in the basis of R and L photon modes. 
As a result, absorption of such a photon leads to valley entanglement of excitons or free carriers. 

Valley entanglement of photoexcited excitons can be verified by measuring the polarization of their PL (Sec. Ill) 
or fluctuations of the photocurrent under an applied in-plane DC bias (Sec. IV). 


II. ENTANGLEMENT OF THE EXCITONIC ENSEMBLE IN A CAVITY 


To visualize the effect of entanglement, here we consider a toy problem keeping only the light-matter interaction and 
neglecting all processes leading to decoherence and losses. In the next section we will include incoherent processes, 
finite bandwidth, fluctuations, and losses. 



FIG. 2: (Color online) Valley entanglement of excitons in K and K’ valleys of a MX 2 monolayer in a cavity geometry. 
Consider a MX 2 monolayer placed into a cavity which supports identical R and L modes of the same photon energy 
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fujj equal to the exciton energy, see Fig. 2. The quantized cavity field is given by 

E = Eij(r)ci^ + Ejj(r)cJ^ + Ei(r)cL + Ej^(r)c^, 
Eii,L(r) = e±Efl,i(r), |Efl(r)| = |Ei(r)| = E(r), 


where the spatial mode distributions are normalized according to 

f s{r)E^{r)(fr = AttHuj, 

Jv 


( 1 ) 

( 2 ) 


where £(r) is the intracavity distribution of the dielectric constant, which may include also the monolayer substrate 
and is assumed to be frequency-independent for simplicity. It can be easily generalized to dispersive media, see 
e.g. [Sllin]. Note that the normalization ([^ takes care of both the electric and magnetic field energy, because for 
stationary fields in a cavity or for periodic boundary conditions one can prove Jy B^{r)d^r = fy e(r)E^(r)d^r. The 
field satisfying Eqs. Q and Q corresponds to the Hamiltonian 

Hf = Hr +Hl = + tiioc^^CL, (3) 

where photon creation and annihilation operators c^,c^ and cr,cl satisfy standard commutation relations [8]. We 
will assume the fields Eir^rIt) and Efl_i(r) to be uniform on the monolayer and omit the position argument. 

We will treat the radiation as quasi-monochromatic and assume that there are a total of N pairs of up and down 
exciton states in a monolayer of area S that can be excited by photons of a given energy. More precisely, TV is a 
number of electron-hole or exciton states within the interband transition linewidth. For a broadband radiation N 
would be a number of states within the radiation bandwidth. 

To simplify the notations, we will denote (f) excitons with latin indices and (J.) excitons with greek indices. The 
Hamiltonian of an ensemble of excitons is given by 


N N 

iJe = ^ E E I’ ( 4 ) 

i=i ?=i 

where the ground state | 0 y 5 )(no exciton) corresponds to zero energy, whereas the state describes an excited 
exciton of energy hoj. The polarization operator of the system can be written as 

N 

Pe = E(dt|l.)(0,l+d;|0,)(l,|) 

N 

where the dipole moments 

dt — d(io)t “ ^-d, d^ = d(io)^ = d(oi)^ = d| = e+d*. 

The Hamiltonian describing interaction between the field and the particles in the electric dipole approximation is 

V = -Pe{ERCR + E*jiC^j^ + ELCL+Elcl). (6) 

Consider Schroedinger’s equation 

ini> = {HR + HL + Hr + Hi + V)'i'. (7) 

We choose the following initial state as a product of the excitonic, 4'e(0), and field, 'l'y(O) wave functions: 

N N 

i=i «=i 

+C'l(0)|1l)|0_r)), 

|C«(0)p + |Ci(0)p = l 


( 8 ) 


4 


in which excitons are not excited and the single-photon field has an arbitrary elliptical polarization. The particular 
case of a linear polarization corresponds to = |C'Lp(see Appendix B). 

The state at an arbitrary moment of time at the same energy huj which is conserved as a result of field-matter 
interaction is 

N N N 

N N N 

+ n n |0,)|0fl)|0,.)e-“‘ 

N N 

N N 

+ Ci(t)iiL)iOfi)nio,)nio?)®”“‘- ( 9 ) 

f=i ?=i 

Substituting into Q, we obtain the equations for coefficients 


thCn = - (d^Ejj) Ef=i C,, ihC, = - (d^Efl) Cr 


ihCL = - (d;E2) Ef=i Q, = - (d;Ei) Cl 
dtE*« 


( 10 ) 


After introducing the complex Rabi frequency 


= n = |r2|e*”, Eqs. (101 can be written as 


iCR = -n c„ ic, = -n*CL 


iCl = -n* Cj, iQ = -ncL 


Taking the second derivative yields 


Cfl,L+iV|flpCfl.L =0. 


( 11 ) 


( 12 ) 


For the initial field in the linearly polarized product state, which is the Bell state (131 in the circularly polarized basis 
(see Appendix B), i.e. when 


CitCO) = ^ 


we obtain the solution 


^ cos(y]V|fl|f) ^ 

Ch =- ^ - ,Cl = ± 


72 

cos (TiVifllt) 


72 


Substituting (141 into (10), we can obtain for 


C,=-i , _ 

^ 7^ 


e”*® ;_ e*® /_ 

sin (7/V|n|f), sin (7iV|fl|f). 


(13) 


(14) 


(15) 


Therefore, the solution of Eq. ([^ when the field is initially in the Bell state lH) and excitons are in the ground 
state has the form 


N 


t=»-*"■ n i“>> n 

i=i «=i 


± |0_r)|1l)) ± fe *‘^‘|0ii;)|0i,) 


72 

sin (7)V|r2|f) 


N 




Vm 




j = l 


N N N N N 

n io.) n io«) ± e*® E 17 ) n 107 n i^-) 1 ■ 

i=i j=i 


( 16 ) 
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For '/N\fl\t = ^ the state (161 corresponds to the Bell-type entangled state of excitons: 




N 


N N 


Eii^)nio*)nio?) 


N N N 

4=1 i=i 


(17) 


This result is intuitively expected: the excitons in each valley couple directly to only R- or L-component of the 
linearly polarized photon state, but these components were entangled with each other, so the entanglement passes 
on to the excitonic system. Clearly this mechanism of entanglement exists only for a quantum incident field. Since 
there is no decoherence or loss, the energy and entanglement oscillate back and forth between the excitons and the 
photon field. If at time t = 0 the cavity contained exactly one photon of a linearly polarized field, then after the 
time t = it will be absorbed by an ensemble of excitons. Since the photon energy is equal to the energy of 


2\^N\n\ 

one exciton, one could (wrongly) assume that this photon creates one (up or down) exciton which will eventually 
recombine, so the reemitted photon will be in either R or L state with equal probability. However, from the exact 
solution (161 we see that the reemitted photon will be linearly polarized. Interestingly, it was predicted in m that 
the spatial structure of a single-photon field should be preserved after its absorption and reemission by an ensemble 
of atoms due to the entanglement of their states. Here we obtain a conceptually similar result for the polarization of 
the field. 

We also note for the subsequent discussion that the state 0 satisfies the following condition for any pair of (f) 
and (4,) excitons: 


— 0, (18) 
where P(ii)j,P(ii),^ are operators of the upper state population: = |lj)(lj|, /5(ii)5 = |l^)(lj|. 


III. ENTANGLEMENT OF THE EXCITONIC ENSEMBLE IN A TRANSMISSION GEOMETRY 


A. The Heisenberg-Langevin formalism 


In this section we consider the response of excitons in a MX 2 monolayer to an illumination by a stationary flux 
of photons, including the effects of relaxation, losses, finite spectral bandwidth, and fluctuations, having in mind a 
generic geometry of Fig. 1. We will solve the Heisenberg-Langevin equation for the density operator of an ensemble 
of (t) or (1) excitons (see e.g. pi IT^'): 


dprnn 

dt 




^mn “f 


(19) 


Here prnn(,'^,t) is the Heisenberg density operator, Rmn the relaxation operator, Fmn is the Langevin noise operator 
describing fluctuations in an excitonic system. The Heisenberg density matrix operator can be determined in a number 
of ways. One can use the projection operator pmn = \n){m\ [HI US]- Another way is to use operators of annihilation 
and creation of quantum states, pmn = al^dm Emil- Both approaches lead to quite similar results. We will call the 
operator matrix pmn a density operator regardless of its representation. The Hamiltonian is 


hmn — bFn^nm — 6)5 


( 20 ) 


where the monolayer is located in the z=0 plane, Wn are eigen values of the unperturbed Hamiltonian of excitons. 
The density operator determined from (19) allows one to calculate the spatial density of any physical quantity A(r, t): 


A{r,t) = {'ifo\AnniPmn{r,t)\'^o), 


where Amn is the matrix element of the operator A corresponding to quantity A and j'l'o) is the initial state function 
in the Heisenberg picture. If the initial wave function is normalized to unity, (4'o|'ko) = 1, the normalization of 
the density operator of an extended system is J2n Pnn{i^,t) = N{r), where N{r) is the density of particles which is 
assumed to be constant in this case. For a monolayer the quantity iV(r) = N{r±) is a surface density, where rj^ is 
the radius-vector in the monolayer plane. 
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The Langevin noise operator satisfies the general conditions = Fmn j {Fmn) = 0 , where (...) means averaging 
over the statistics of a noise reservoir. Other properties of the noise operator are determined by the properties of the 
relaxation operator. We will define the latter in the simplest form through a constant relaxation rate: 


Fmn — ^mnPmn: TTl ^ 71 

Fmm — ^ninPnn 


( 21 ) 


The definition (21) implies that the relaxation operator is delta-correlated in time. Using generalized Einstein relations 


and assuming also delta-correlation in space, one can derive the following relationship for the noise correlators 
for a two-level system nans]: 


{F^n{t, r±)Fmn{t'= {2'ymnPmm + Rmm)S{t - t')6{r± - 

{Fmn{t, r±)F^^{t'= {2'ymnPnn + Rnn)5{t - t')S{r± - r'j_). (22) 


In a two-level system, we denote the transverse relaxation rate as ■jmn = 7 io = 7 - Assuming that the exciton 
system is not too far from the equilibrium and the Fermi level is far below the exciton band edge (so that the number 
of excitons is always small), we can also obtain 


Rii — —Roo ~ —Tpii. 


(23) 


Under these conditions and for 7 T we can also take Rmm 


0 and Rnn « 0 in Eqs. (22) and (44) as in [3. 


B. Response of the MX2 monolayer to an incident field 

Now we derive the optical-frequency current induced in the valleys and the properties of the reemitted field. We 
start from the following ansatz for the field operator in the monolayer E(z = 0) = E(0): 


E(0) = 


E 


(-) 


d+) 


(-) 


; E = e+EnF Q-El, E = 


(-'")’■ 


(24) 


Next we will use the equations for the density operator (19) and (20). They can be split into two pairs of equations: 
for t FR and I +L states. These pairs can be coupled via an entangled initial state as we showed in the previous 
section; however, this does not affect the form of the equations. 

Introducing slowly varying amplitudes of the off-diagonal density operators (’’coherences”). 


P(io)t7 — o'(io)t7e 


we obtain the density operator for an R-and L-polarized surface current: 

hi = (e±jT7e—‘ + e4,,e“*) 

= (e±J7^i3-(io)t,;e -h e;pj^7(7(io)^7e“‘). 


(25) 


(26) 


where we introduced the notation j(io)t = Jt = j(io)^ = ji = iujd*. 

Assuming that the spectral bandwidth of the radiation does not exceed the relaxation rate 7 we obtain the stationary 
solution of Eqs. (19), (25), (26): 

» iud* ~ 

Jl — j-_, (.P(00)t ~ P{ll)l)^R - 


H'y 


7 

iujd . 


Ji - -^{pmi - P(iiu)Fl 


/5(ii)t - (jIFr + E^Jt) 

P{ii)l Yhujj 


-h 


E, 


(ii)T 

r ’ 

E(ii); 


(27) 


where = E(xq)^^^. We will assume that fluctuations in different valleys are not correlated: 


{F\F^) = (E(ii)^F(ii)^) = 0. 


(28) 
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For a single photon flux we can safely assume that the surface density of excited excitons is much lower than the 
maximum density of exciton states N determined by the spectral bandwidth Aw and the density of states. Then, 
taking into account the normalization condition 


P(oo)t + P(ii)T — P(oo)^ + P{ii)i — ^) (29) 

we can replace P(oo)^ — P(ii)^ —>■ N and P(oo)^ ~ A in Eqs. (27), which gives 

4717 7 4717 7 

where = -U- jg the surface cooperative frequency squared which has the dimension of [cm/sec^]. 

Equations (30) and the last two Eqs. (27) give a complete description of the response of electron-hole or exciton 
states in a MX 2 monolayer to the field (24). In particular, the polarization of the surface optical current excited in 
the monolayer and of the radiation emitted by the current can be found by calculating the Stokes parameters: 


— (jljy}, 


(31) 


where 


and the brackets (...) mean averaging over the statistics of the reservoir and over the quantum state. 


C. Self-consistent optical field in the monolayer 


Next, we need to relate the unknown field Eq. (24) in a monolayer to an incident field. For simplicity we will 


consider a ID propagation problem, corresponding for example to a wide enough beam. If the monolayer is located 
in the z = 0 plane on a substrate with dielectric constant e = v? for z < 0, the electromagnetic fields incident on the 
monolayer from both directions can be written as 


(i) Pump: 


E,(z<0) = E,(e+e^’^^-^’^*CR + e_e 


+ e_e 


ikz — iujt 


Cl + e+e 




(32) 


(ii) Vacuum noise: 


E,(z>0) = F;o(e+e-*'="-“‘cyfl + e_e 


ikz-\-iujt 


-VR 


e_e 


—ikz — iujt 


cvL + e+e 


ikz-\-iLijt 


-VL, 


(33) 


Here = y/27rliw/e and Eq = ^/2TThui are the normalization amplitudes. In a more general and realistic case, the 
incident field is not exactly monochromatic but occupies a narrow frequency band Aw ^ w. In this case the creation 
and annihilation operators are defined for slowly varying amplitudes of the Heisenberg operators (see also [51 110)1. 
Note that the operators cr and cvr^ and other similar combinations (same for cl,cvl) correspond to the field modes 
with different wave numbers: k = \k\ and k = — \k\. Therefore they always commute. 

Normalization amplitudes in Eqs. (32 33) correspond to a unit quantization volume V = 1, i.e. the dyadics c^rCr 
and c^cl are the operators of the photon density (see [HlinKlT]), and the quantities like (c^rCr) = ('I'o|Cfl;Cfl|4'o) and 
(crCl) = (4('o|cJ^cl|'I'o) have the dimension of [cm“^], where 'kg is the initial state function present in the Heisenberg 
picture, which is normalized as (4'o|4('o) = 1 . 

When calculating relative average values, it is sufficient to take the field as monochromatic and use the explicit initial 
state 4 / 0 . For example, if the initial field is in the linearly polarized single photon state equivalent to a single-photon 
Bell state in the circularly polarized basis, 'kg = 'ki. 


= 


|li?)|0L) ± |0ij)|lL) 

V2 


(34) 


















the averages have the following properties: 


(CrCr) = (CiCi) = ±(4 cl), (4cl) = (4cfi); 

{CrCRC^L^l) = 0. 


If the incident field corresponds to the other Bell state, e.g. 

|0i^)|0L) ± |lfl)|lL) 


= 


V2 


(35) 


(36) 


or if it is a state with the same average energy Huj which is not entangled in the circularly polarized basis, e.g. 


we obtain instead of Eqs. (351 that 


'I'd 


\0r) + e^'^llR) |0L) + e''^|lL) 

V2 "" 


(crCr) = {c{cl), {c’rCl) = = 0; 

(CrCRcIcl) = (4ci?)(4cL)- 


(37) 


(38) 


Whenever it is important to include Langevin noise terms, one has to take into account the finite spectral bandwidth 
of radiation Aw ^ w. In this case the properties of commutators and correlators are determined taking into account 
the density of photon states. For example, for a paraxial beam with the aperture cross section S± we obtain following 
pun] that 


27rcS' 


-6{oj — oj'). 


(39) 


where = 4 ^ Co-^c 


*dw, n is the refractive index, and the subscript a denotes the polarization of a normal mode 


of the field. For correlators of spectral components of the field satisfying (E) = 0 we have 

-S'-L-(cLw) = 

n Ztt 


(40) 


where the dimensionless quantity Na-c^/2Tr determines the power Pa-u, = hui{N„^/2 t:) incident on the area Sa_ per 
unit frequency interval per unit time for a given polarization, so has a dimension of energy; for the vacuum field 
= 0. Here we won’t present cumbersome expressions for the state functions of multimode fields which are reduced 


to states (341, (36), or (37). We will only give some correlators needed for the subsequent derivation. In particular, 
conditions (35) and ([^ correspond to 


(CrCr) = (4 Cl) => Nr^^ = Nlui 

(CrCl) = {c^rCr) = ±(4cd) ^ {c'n^CLuj') = {Cr^CRuji) = ±{c'h^Cr^>) 
{c^rCl) = (4ci?) = 0 => = {c\uj^Rui') = 0- 


(41) 


Note that since a duration of a single-photon ’’pulse” is of the order of the maximum power of a single 

photon illumination regime within the spectral bandwidth Aw of of the order of SwAw ; for a higher power individual 
photons overlap. Here we don’t consider the possibility of using ’’multiphoton” fluxes of quantum correlated fields. 

After performing a standard calculation of the operators of the reflected and transmitted fields, we obtain the 
following expression for the field operators Er^l used in Eq. (24): 


Er^l = 


2Eo 47rj;,^ 

{cr,l + cvr,vl) - 


n -I- 1 


(n -I- l)c' 


(42) 


The last term on the right-hand side of Eq. (|42|) gives rise to the the collective superradiant relaxation of the excitonic 
ensemble with the rate 0.sr = 




c L 


(n -I- l)c 


20]. Taking into account Eqs. (30) it is easy to find that this term can be 


neglected provided 7 ^SR- In this case one can neglect radiative corrections (back reaction) when calculating the 
current excited in the monolayer. 
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D. The effect of radiative corrections 

To study the effect of radiative corrections on the excited state populations (pn) it is convenient to introduce a 
space-time spectrum of noise, 


p — 
m.n. — 


q,v,;mn 


iqrj_-TOt T 72 p _ pt 

c LLULL X q^v^rnn — ^ —q —^ 


(43) 


we can obtain from (22) the correlators for the spectral components of the noise operator: 


{^q ,v,mn^q' ,v' ^rnn) — ^^3 (^mn{Pmm) + V )( 5 (q Q ) 

{Rq,v,mnFq/^yi-^rnn) ~ \ lmn{Pnn) + '^{Rnn) \ — V )l 5 (q — q ). 


(44) 


Note that the harmonics of the current with wave vectors q > to jc cannot excite propagating modes and therefore 
cannot give rise to radiative losses. Therefore we can put q, q' < lo/ c Eq. (44) and, using Eqs. ( |27|30|42|44D , after 
lengthy but straightforward calculations, obtain an estimate of the radiative relaxation rate of populations induced 
by noise: 


Vl 


rad 


2w|dp 

ch 



(45) 


As expected, this rate corresponds to the inverse spontaneous emission time. 

To summarize, when calculating the excited state populations {pii) the radiative corrections can be neglected if 
2nsR, ^rad E. Taking into account the back reaction effects will make the derivation more complicated, and the 
result will amount to renormalization of the relaxation rates in the final expressions: 


7 — > 7 -|- flsR, r —>■ r -|- 2^1sr + ^rad- 

To keep the derivation more streamlined, we will ignore the radiative corrections and assume that 

2£’o ^ 

Rr,L = - —\Crl + CVR,VL)- 

n + 1 


(46) 


E. The polarization of the photoluminescence from excitons excited by a single-photon field 

The initial (Heisenberg) state of our system can be represented as 

l^o) = l^n)l^y)l^(i?-i-L)), l^'n) = |0t)|0 4.), 

| 4 'y) = | 0 yfl)| 0 vL). ( 47 ) 

First let us neglect the contribution of the Langevin sources. Consider an incident field in the Bell state 


14- 


R+L 


) = 


(48) 


which corresponds to a single-photon linearly polarized field |l2,)|0y) or |0a;)|ly) respectively (see Appendix B). If we 
substitute expression (46|) in Eqs. ([30|) and ([sT]), then after performing the averaging in Eqs. ([sT]) all dyadics that 
include the vacuum fie 


CrCr, 


c{cL, 


C^rCL, 


d operators will become zero. A non-zero contribution comes only from averaging the dyadics 


and Cj^cr: 


(«'ol4cfll^o) = (4 'o|4cl|«'o) = 

(4'o|4cl|^o) = («'ol4cfil^'o) = ±^. (49) 

As a result we obtain ^ 0, Sy = 0 for = 4'+ and Sy ^ 0,3^ = 0 for In other words, the 

reemitted photons will have linear polarization as in the cavity geometry described in Sec. B. At the same time, the 
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reemitted photons will be unpolarized, Sx = Sy, if the incident field corresponds to other Bell states of a single-photon 


energy, Eq. (36), or if it is a single-photon state which is not entangled, Eq. (37) 


So what is the role of the exciton entanglement here? The fact that a linearly polarized radiation remains linearly 
polarized after interacting with a system of equal numbers of ’’left” and ’’right” rotators (such as K and K’ valleys in 
M 0 S 2 ) does not by itself constitute the evidence that left and right subsystems are entangled. Indeed, the same result 
can be obtained for a classical field [ 6 ] which cannot entangle the non-interacting quantum systems; see Appendix A. 
However, for a quantum incident field, e.g. for a single-photon field the conservation of the linear polarization in the 
reemitted radiation is possible only if left and right excitons are entangled. Therefore, in this case the degree of the 
linear polarization of the reemitted field can be used to verify the exciton valley entanglement. 

In order to get an explicit expression for the noise-induced depolarization of the reemitted field, we need to sum 


over all modes in a spectral bandwidth Auj. If we take the Langevin noise terms into account in Eqs. (30) and use 
Eqs. (39) and (40), the Stokes parameters for the initial state = 4'+ are 




\27r7(n -|- I) 


-b Nluj)Au 


2'!tcS± 


+ A, Sy — A, 


A=‘^{PIP, + PIP,). 


(50) 

(51) 


When calculating the value of A we take into account only the spatial harmonics of the Langevin noise within the 
spectral bandwidth Aw that correspond to a paraxial beam of the aperture S±. We will also assume that the relaxation 
of coherence is much faster that the recombination rate, 7 3 > T, which is typically the case in semiconductors. Using 
also Eq. (441, we arrive at 


A = 


(4>o|p(n)^|4'o) -b (4>o|p(n);|4'o) 
2Tr‘jS± 


Sx a -b 1 


(52) 


a = 


2 iV(iVfl„ + Nl^) 


C7 


'^c± 


(n -b 1 )^ 
2n 


In order to find the fraction of the unpolarized field a we need to calculate the operators P{ii)^ and P(ii)\. in 
Eqs. (1^. Substituting (ISOl) and (|46|) into the last two Eqs. (p^ it is easy to find that the resulting noise-dependent 


terms are linear with respect to the Langevin noise operators, so that they vanish after averaging over the reservoir. 
As a result, we obtain 


/3(ii)t = X(4 + c\^r){cr + cvr), 
Piii)l = X{cl -b c|)^)(cl -b cvl), 


(53) 


where x = 


^c L 


\Eo? 


Trjr(n + l)^hui 


After summing over the modes within the bandwidth Aw and taking into account 


Eqs. (39),(40), we get 


('kol^Dtl'f'o) = ('kolP(iiul'ko) = 


27r 


Substituting Eq. (54) into (52) yields 


Aw I 
^ 2NS_l ’ 


(54) 


(55) 


Where the ratio N/Auj is determined by the exciton or electron-hole density of states. Equation (55) has a simple 


physical interpretation: a fraction of the unpolarized field in the reemitted radiation is equal to the ratio of a number 
of photons incident on a monolayer during the lifetime of an exciton to the total number of exciton states within 
Aw. As is clear from Eq. (55), a scales as one divided by the number of exciton or electron-hole states within the 
bandwidth T ~ 10® — 10^® s^^ so assuming a standard 2D density of states for parabolic bands and the beam 

aperture S± > ^ 10“® cm“®, a is smaller than 10“®. 
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F. Correlation properties of photoexcited excitons 


It is easy to show that in the case of the initial Bell state of the field given by Eq. (34) the correlator for the density 
operators (531 obeys the following property: 


('^'o|P(ii)tP(ii)tl'I'o) — 0, 


(56) 


which coincides with the corresponding relationship ([^) for the solution to the Schroedinger equation for a cavity 


field. It is important that when Eqs. (28) are satisfied, t 


le condition (56) is not affected by the Langevin noise. 


For unpolarized but still single-photon quantum fields (|36[),(|37[) we obtain instead of ( |56[ ) that 

('f'o|P(ii)tP(ii)tl'^'o) = ('f'olP(ii)tl'^'o) X (^o|P(ii)^l'f'o)i (57) 

which is the same as for classical fields. Therefore, Eq. ( [5^ is a unique property of the valley-entangled exciton system 
created by an incident quantum field in the initial state (48). Excitons created by a classical field or by quantum 
fields in the initial states (36|37) do not satisfy Eq. (56). 


IV. PHOTOCURRENT IN A VALLEY-ENTANGLED ELECTRON-HOLE SYSTEM 

Quantum correlations between photoexcited carriers should manifest themselves in the fluctuations of the photocur¬ 
rent or photovoltage under the applied bias. A DC electric field applied along the monolayer can lead to separation 
of photoexcited electrons from holes while not affecting their valley index and any possible entanglement with respect 
to the valley degree of freedom; see Fig. 3. This will create a photocurrent or photovoltage depending on the way the 
detector is wired in an external circuit. Note that breaking the binding energy of excitons may require a very high 
electric field. In experiments [7] the charge separation most likely originated from the metal-semiconductor contact 
regions with a high built-in electric field. 



FIG. 3: (Color online) A sketch of photocurrent detection. Under illumination with a linearly polarized single-photon field 
and in the presence of an in-plane DC electric field Edc there is charge separation in the direction along Edc which gives 
rise to a photocurrent Jp if a monolayer is contacted and included in a proper external circuit. There is no photocurrent or 
photovoltage in the transverse direction. 


A lot of attention has been recently devoted to the valley Hall effect which generates photovoltage in the transverse 
direction under a circularly polarized excitation El El- However, for a linearly polarized excitation there will be zero 
Hall voltage with equal amounts of positive and negative charges accumulated on the sides as sketched in Fig. 3. 
Here we consider the longitudinal photocurrent along the direction of the applied DC field. Specifically, the current 
detector will measure the quantity Jp = (d^ol jp|4'o), where the averaging is taken over the initial state (47) and Jp 
is the photocurrent operator: 


Jp — ’7(P(ii)t + P(ii)4)) 


(58) 


where ry is a certain coefficient. Obviously, equal amounts of photogenerated electrons and holes from K and K’ valleys 
will contribute to the signal. When characterizing quantum correlation properties of carriers, the quantity of interest 
is not the current itself but current fluctuations defined as 


5Jp = ^(J|,)-(Jp)2. 


(59) 
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Using the initial quantum state defined by Eq. (471, we compare the current fluctuations due to the following three 


states of an incident field that have the same average energy: 

(i) Unpolarized field: 


|0fl) + e*1li?) ,, |0i)+e*’^|U) 
'^R+L = 'i’D = -^- X 


(ii) Circularly polarized field: 


V2 


'^R+L = '^R= |lfi)|0L). 


V2 


(60) 


(61) 


B): 


(iii) Linearly polarized field that can be represented as an entangled state of two circular polarizations (see Appendix 


4'i^+L = 4'+ = 


|lfl)|0L) + |0fl)|lL) 

72 


(62) 


To calculate the photocurrent fluctuations, we need to evaluate the averages of both the excited-state density 
operators in Eqs. (53) and their products. The first of them gives 


P(ii)t — xi^^R^R T ^vr^r “b ^ffivR T ^vr^vr)] 


At 




At 


(63) 


P(ll)tP(ll)T — X -|- C^trCvrCrCr 

~^C^rc\iCvrCvR -f CyRCyj^CyR Cyj^CRC^j^CR 

+CRcl^CRCyR + cl^CRCyj^CyR + cl^CRCRCyj^ 

+CyRCyj^C^j^CR + Cyj^Cyj^CyRCR -|- 
+CyRCvRCRcl^ + Cyj^CyRCyRC^j^ + C^RCRCyj^CyR 

+CyRCyRCVRCR + Cy RCy j^CyRC^f/). 


(64) 


The same result can be obtained for the population of the states after replacing f with I and R with L. After 
averaging (63) and (64) over the initial quantum state Eq. (47) and taking into account the commutation relations 
we obtain: 


= X{^R+L\c]iCR\^R+L), {p(ii)l) = x{'^R+l\ClCl\'^R+l), 

{P{ii)tP{iin) = X^R+lIc^rCrI'^r+l) (l + {^VR\cyRc\,p.\QyR)^ 
(P(ll)t/5(11);) = {'^ R+l\c\cl\'^ r+l) (|l + (0yL|cvLc[.^|0yL)^ . 


(65) 


( 66 ) 


As is clear from Eqs. ( 661 , vacuum fluctuations of the field amplify the fluctuations of the photocurrent. In the absence 


of an incident field Jp = 6Jp = 0 as expected. It also follows from from Eqs. (65), (66) that the following relations 
are true for all three states of the field: 


(P(ii)t) — (P(ii)4) — Xj (P(ii)tP(ii)T) — (P(ii)t/5(ii)t) — 2x^- 
Next, we calculate the averages for ’’mixed” products P(ii)^P(ii)i and /5(ii)^/5(ii)^: 

(P(ll)tP(ll)4) = {P{ll)lP{lin) = ^X^ R+L\cRCRC'^CL\'i’ r+l) ■ 


(67) 


( 68 ) 


Proceeding in the same way as above, for an unpolarized field ( [6()| we obtain (/5(ii)^P(ii)^) = (p( 
whereas for polarized fields (61) and (62) (/5(ii)^/5(ii)^) = (/5(ii)^P(ii)^) = 0. This result is trivial 
polarized field (61) since the Tarter does not interact with excitons. Eor a linearly polarized field 
is not obvious and is entirely due to the entanglement of and ”t” excitons. 

Using Eqs. (67) and (68) to calculate the current (|58| and its fluctuations (59) we obtain SJp 


ii)t) ^ (P(ii)t)’ 
for a circularly 
(1^ this result 


illumination with an unpolarized field and SJp = Jp in the case of a linearly polarized incident field, 
of fluctuations by 72 is the direct consequence of the interference within the valley-entangled state 
carriers. 

In conclusion, we have shown that excitons in K and K’ valleys transition metal dichalcogenide 
be efficiently entangled by interacting with linearly polarized single photons. Valley entanglement 
polarization of reemitted photons and squeezing of the photocurrent fluctuations. 


= 72 Jp for an 
This reduction 
of photoexcited 

monolayers can 
leads to linear 
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V. APPENDIX 

A. Conditions for entanglement with classical fields 

Consider two quantum systems described by generalized coordinates qi and q 2 . Their wave function 'i'{qi,q 2 ,t) 
obeys Schroedinger’s equation 


ih'ii = H'^. 


(69) 


Assume that these two systems are not coupled with each other directly but inte ract with a classical electromagnetic 
field described by a classical variable u{t). In this case the Hamiltonian in Eq. (69) can be written as H = Hi {qi, u{t)) + 
H 2 {q 2 ,u{t)) where the operators Hi and H 2 act only on the functions of variables qi and q 2 , respectively. Then, 
substituting 'k = il;i{qi,t)'ip 2 {q 2 ,t) into Eq. (69), we obtain 


1p2 (^iflTpl — -I- '01 — 77202^ = 0. 


(70) 


If the quantum systems were not entangled at the initial moment of time t = 0, i.e. 'I^(qi,q 2 , 0) = 0i(gi, 0) x 02( 92 ,0), 
Eq. (70) splits into two independent equations for each system: 


ihipi = Hitpi, ih'tp2 = 77i02, 


(71) 


and the solution will remain in the form of the direct product 'i>{qi,q 2 ^t) = ipi(qi,t) x 02 ( 92 ,i), corresponding to 
unentangled systems. 

Entanglement may appear if the classical field gives rise to the interaction Hamiltonian V(u(t),qi,g 2 ) that directly 
couples the two systems, for example: 

=-ffi(9i) + ^ 2 ( 92 ) + X'«(09i92, (72) 


where y is a coupling constant. This particular example corresponds to the Hamiltonian describing generation of 
entangled photons in a medium with a second order nonlinearity as a result of parametric frequency conversion [S]. 


B. Relationship between XY and RL basis states for photons 
Consider the quantum field in vacuum in a quantization volume V: 

E tt' ikz — iujt^ I —ikz-\-icjt'■i . ,, ikz—iut'' . ,, —ikz-\-iujt 

= ho (xoe Cx + xoe ^ + yoe Cy -I- yoe ^ 

where Xo,yo are unit vectors along x, y coordinate axes, 

/ sin kz d^r = / cos kz d^r = / sin 2kz d^r = / cos 2kz d^r = 0, 
Jv Jv Jv Jv 


yJ ’ 


(73) 


(74) 


En — 


2'Khu} 


V 


Let’s expand the field given by Eq. (73) in terms of circularly polarized orthogonal modes: 


E = Eo [e+e 


ikz—itot 


cr + e_e 


— ikz-\-iujt 


e_e 


ikz—iujt 


Cl + e+e 


— ikz-\-iLjt 




(75) 
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where e± = ■ Comparing Eqs. (731 and (75) we obtain 


_ cr± Cl „ _ c^TCy 

^x,y — 5 ^R,L — 


(76) 


V2 ’ V2 ■ 

Next, we introduce the vacuum state |0 e); (OeUOs) = 1. By definition of creation and annihilation operators, 

Nr / . \ Nr 

b»J (4) fe) 


lOv 

|iV,,)|iV,) = |iV^)|iV,.) = 


Expanding the linearly polarized photon state 


V^rINlI 


«' = IU|0,) 


(78) 


in the basis of circularly polarized modes, we obtain: 


X! ^Nr,NlWr)\Nl)- 

N ji—0 ,Nl—0 


Taking into account Eq. (76), one can obtain from Eqs. (78) and (79) that 


(ct 

^Nr,Nr = (N‘i^|(fVL||l2;)|0y) = (0s|— ^^===^4|0e) 


= (0j 


{cr)^^{cl)^^ (4 +Cl) 




|0e) = 


0 


if fVfl + iVL ^ 1 


IIrOi, 


= ^OhIi, = ^ if + Nr = 1- 


(79) 


(80) 


As is clear from Eq. (80), the product state |l2;)|0j,) in the basis of linearly polarized modes is one of the Bell states 
in the basis of circularly polarized modes: 


= 


|li?)|0L) + |0i?)|lL) 

■ "V2 


(81) 


The relationship |l3;)|0y) = can be also obtained in a less formal way from average values of squares of cartesian 
components of the field. Indeed, for the field given by Eqs. (73), (74), and (75) one can obtain the following expressions 
for the averages: 

^ =2hw, t 44 d^r\ = Huj (1 + l] , ^ ( f EyEy d^r) = (82) 


47r 


Using the relationships 


47r 


2y ’ 4^ \Jv ^ ^ 


Ex = —F= 1 e 

V2 V 


P / 

~^0 I ikz—iujt ^ j_ ^ikz—icot' 


CR + e 


ck + c CL + e 


—ikz-\-iujt 


). 


(83) 


4 - ^ (*e* 


kz — iujt' 


— ikz-\-iojt ■ ikz—iujt' 


-cr - le -Cr - le-- -'‘Cl + ie ^ 

we can express the averages through creation and annihilation operators of circularly polarized modes: 


(84) 


(85) 


Huj 


ExEx d'^r) = '— (cfi4 + CflCfi + clC^ + c 4 l + “^crc^ + ^ErCl) , 


t A , 


U , od 


( 86 ) 


EyEy d-^r) = ^ (c_r 4 + CflCfl + cz,c^ + “ 2 (crc^^ + 24cl) ) 


(87) 
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The first of Eqs. (82) can be satisfied by any states of the type 

^ = A\1r)\0l) + B\0r)\1l), 
if 

\A\^ = \B\^ = i 


( 88 ) 

(89) 


Taking into account the last two equations in Eq. (82), we obtain 


{\A\^ + \B\^) + ^{AB*+A*B) = ^, 


{\Af + \Bf)-l{AB*+A*B)=^-, 


which yields 


A = B = 


71 ’ 


which is equivalent to the result in Eq. (80) up to insignificant common phase (j). 

It is easy to see that if we start from another linearly polarized product state |0a;)|ly) instead of Eq. (78), we arrive 
at another Bell state: 


72 


(90) 
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